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Abstract 
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1. Introduction 

Heavy Quark Effective Theory (HQET) is an effec¬ 
tive theory for QCD, the theory of strong interactions, in 
the limit where quark masses are large and other scales, 
such as momenta are kept fixed. Understanding this 
limit is of great interest per se. In addition a control of 
HQET is very useful to arrive at phenomenological pre¬ 
dictions for B-meson properties and qualitatively also 
for D-mesons. In particular B-meson decays need to be 
understood better in order to further constrain the flavor 
sector of the standard model of particle physics. 
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In this article we give an introduction to HQET 
with an emphasis on its full non-perturbative formula¬ 
tion. Mostly we remain with the general ideas and an 
overview of the present status and results. For more de¬ 
tails concerning the basics as well as the phenomenol¬ 
ogy we refer to the literature JT| [2, 3] [4j [5) with the 
non-perturbative aspects and particularly the discretisa¬ 
tion on a lattice covered in the last reference. Here, on 
the other hand, we give a more complete discussion of 
the status and of non-perturbative computations and the 
challenges for the future. 

HQET, as discussed here, is an effective field theory 
for the low energy physics of energy levels or transition 
matrix elements with a single heavy quark or anti-quark 
in initial and/or final state. We will label these hadronic 
states by H and the quark by h. The latter has a mass 
m h. [[] Usually we think of H as a B-meson, but it can 
be e.g. a baryon with beauty quantum number of one. 
We always consider a rest-frame where all spatial mo¬ 
menta p, are small. The effective theory then yields the 
expansion of observables of QCD in powers 

(IP/I /m h ) n , (A /m h )", 0 mj/m b y (1) 

where h is a heavy quark, while the masses nij of the 
other quarks are considered small. 0 

1.1. Mass scaling and phenomenology 

The HQET expansion is of a theoretical interest be¬ 
cause it describes the asymptotics of QCD as mu —> °o. 
For example we obtain statements such as 

W(qcd({P,■),{'«;)’'» h) = A4 slat ({p 1 },{« ; },m h ) 

x (1 +0(1/m h )) (2) 

with 0(1 /mu) summarizing the terms of eq. 0 The in¬ 
trinsic scale A may be taken to be any low energy QCD 
scale. 

The important content of eq. (|2]) is that it gives the 
large mass scaling of observables AIqcd with mu, in the 
form of the static (lowest order) of HQET prediction, 
usually 

M stat ~ (muf (3) 

with the (not necessarily integer) power s determined 
by counting dimensions and adding anomalous ones in 


1 We use the symbol mu generically when the precise definition of 
the renormalized mass does not matter. 

2 At this stage we are interested mostly in the theory and less the 
phenomenology, where one may ask whether the charm quark is to 

be treated as a light quark or a heavy one. Top quarks on the other 
hand are not considered at all. They are heavy enough to safely be 
considered as decoupled. 


the static effective theory. Understanding this scaling is 
clearly a very relevant part of understanding QCD. 

A second important motivation for studying (and 
computing in) HQET is that the b-quark mass, say in 
the MS scheme at 4GeV renormalization scale, is of 
order 4GeV. It is an order of magnitude larger than 
the intrinsic QCD scale of around A * 400 MeV. In¬ 
deed, if one wants the 0(1 /mO to give a first estimate of 
the numerical size of the corrections (without additional 
pre-factors) this value for A is appropriate in eq. 0 as 
we will see in section [5] Consequently, as long as we 
keep momenta small, static predictions are expected to 
be good at the 10% level and one has an accuracy at 
the 1% level when 1 /mu corrections are included. Thus 
HQET is a very interesting phenomenological tool. 

1.2. Heavy quarks in lattice QCD 

One may still wonder why it is of interest in the con¬ 
text of lattice QCD. The reason is simply that a numer¬ 
ical lattice QCD computation necessarily is done with 
an infrared cutoff l/L through the linear extent L of the 
simulated T x L 3 world on top of the ultraviolet cutoff 
1/a introduced by the lattice spacing a. The accessi¬ 
ble physical energies E, have to be removed from these 
scales, 

1 /L <K Et « 1 /a , (4) 

otherwise properties of the associated states are dis¬ 
torted. 

In table Q] we list the most relevant effects that are at 
the origin of these bounds as well as the errors which 
result from violating them. Since the finite volume ef¬ 
fects are exponential in L, the bound of m„L = 4 is 
rather sharp. However it depends on the pion mass m„ 
which one has in the simulation. Reasonable values of 
m n = 300 MeV ... 150 MeV lead to L > 2.5 fm... 5 fm. 

In contrast discretisation errors only disappear like 
a 2 . Further they vary a lot depending on the quantity 
and discretised action. They simply have to be stud¬ 
ied by changing a and the difficult question is where 
the asymptotic a 1 behavior sets in. From then on a fac¬ 
tor 2 variation in a 2 (or better more) is acceptable. We 
have included our rough estimate where a 2 scaling sets 
in. Together with the required L > 2.5 fm... 5 fm, this 
shows that L/a has to be prohibitively large when the 
b-quark is simulated as a relativistic quarljj while for 


3 We should mention that not everybody in the field agrees with this 
statement. There are lattice QCD computations with quark masses 
very close to the physical b-quark mass and amu < 1. Discretisa¬ 
tion errors are fitted with polynomials in the lattice spacing and these 
representations of the data are used to extrapolate the results to the 
continuum and the physical mass. As an example we cite ED 
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source of errors 

cases 

asymptotics 

% effects 

References 

finite volume effect due to 
particle exchange around 
the periodic space 


0(exp(-m gap L)) 

m n L =3 4 

000 

discretisation errors 
(O(fl)-improved) 

lattice QCD 
with b 

0((aEj) 2 ), 

Ei ~ m n ... m B 

extrapolation with 
a < 0.025 fm 

smalm 


lattice HQET 
without c 

0((aEj) 2 ), 

Ei ~m n ... |p/|, A 

extrapolation with 
a < 0.1 fm 

0EEMD 


Table 1: Effects in lattice QCD computations due to infrared and ultraviolet cutoffs. The lowest particle mass is denoted by ra gap ; in QCD with 
light quarks this is m gap = m n . In the column titled “% effects” we list the condition needed to have systematic errors around the 1% level. 


HQET lattices of size L/a = 32 ... 64 seem sufficient. 
For this reason lattice HQET is a very attractive phe¬ 
nomenological tool. 

1.3. Defining effective field theories beyond perturba¬ 
tion theory 

We now describe the general concept and formula¬ 
tion of an effective field theory. The special features of 
HQET will be mentioned in the following subsection. 
We consider processes in a fundamental theory (QCD 
or the standard model of particle physics - the impor¬ 
tant feature is the renormalizability of the theory) at low 
energy. In particular we first focus on processes (scatter¬ 
ing, decay) of particles with masses of this low energy 
or below it (in HQET also the large mass particles are 
involved as will be discussed soon). In this situation, 
vacuum fluctuations involving much heavier particles 
are suppressed and a true creation of the heavier par¬ 
ticles is energetically forbidden. One therefore expects 
to be able to describe the physics of these low energy 
processes by an effective field theory containing only 
the fields of the light particles m. The leading order 
Lagrangian of the theory is formed first from the free 
field theory Lagrangians and all the renormalizable in¬ 
teractions. Assuming the usual power counting, all local 
composite fields with mass dimension smaller or equal 
to four are allowed. Let us denote the Lagrangian by 
«i? L0 and the Euclidean action is .S' 1 0 = f d 4 xJf LO (x). 
Correlation functions are then defined by the standard 
path integral 

(O)lo = ^~ f e- SL °0. (5) 

Z LO J fields 

with (1 )lo = 1 and O some multilocal product of fields 
such as O = <t>(jc)0(v). In this way we start with a renor¬ 
malizable theory. For a lattice formulation this means 
that the continuum limit of the theory exists when a fi¬ 
nite number of renormalized parameters are kept fixed. 


The continuum limit is then also expected to be univer¬ 
sal, i.e. independent of the specific discretisation. 

Higher order terms in the expansion of physical am¬ 
plitudes (or correlation functions) in 1 / m h are given by 
including fields with higher mass dimension, which is 
compensated by the appropriate factor of the large mass 
in the denominator, 

jf NL ° = yW-, 0 j i = —cb i (6) 

where the parameters oj, are dimensionless. The fields 
contained in the (multi-local) O are expanded in the 
same way as the action, <9 e ff = Oi.o + Omo + . .. . We 
now have to deal with interactions in eq. which are 
not renormalizable (by power counting). However, we 
are only interested in the expansion <t> = ffij, | ( | ) + <lQj'° + 
... of observables <E> in m h l . It is therefore sufficient to 
define the theory with the weight in the path integral ex¬ 
panded, e s —> e~ s {1-5 NLO + ...}. At NLO accuracy 
the expansion is then given by 

d)LO = <0 LO )lo (7) 

onlo = < 0 nlo )lo (8) 

- (<<9 lo S nlo > lo - <D lo )lo <5 nlo >lo) 

and S NL0 = J d 4 x _i? NL0 (x). The term <lA f| L0 (but not 
the individual terms on the r.h.s. of eq. (|8]>) is renormal¬ 
izable with a finite number of counter terms which are 
equivalent to renormalizing the parameters w, (includ¬ 
ing the LO ones). Also parameters in the fields O are 
part of the list of co,. The reader may worry about di¬ 
vergences in the form of contact terms between Jz? NL0 
and 0 LO in (0 LO S m '°)lo- These, however, can all be 
absorbed into the oj,, see D3 and 0. 

Renormalizability is particularly important for a non- 
perturbative evaluation of the path integral in a lattice 
formulation. The continuum limit of an effective theory 
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only exists when we treat the higher dimensional inter¬ 
actions as insertions in correlation functions in the form 
of eq. <§. 

1.4. HQET 

HQET reaches somewhat beyond the situation dis¬ 
cussed above. The difference is that we are interested in 
processes which do involve the heavy quark h at small 
momenta (remember we choose the rest frame prop¬ 
erly). It is therefore not immediately clear, what are the 
degrees of freedom to be kept at low energy: what is the 
complete basis of low energy fields? The answer to this 
question was found in various ways. We here sketch one 
line of reasoning. 

For smooth fields, the Dirac Lagrangian 

«5f = i//(m h + D^y^ijj (9) 


can be split order by order in 1/mh into decoupled up¬ 
per and lower compenent quark field contributions, cor¬ 
responding to the particle and the anti-particle field: 


^ = 

J? h stat + ^(D 

(10) 

+ 

+ 2Sf, (1) + OC-4) 
h h ml 

(11) 

cgs tat _ 

+ D 0 )i//h, 

(12) 

c/) stat 

<Ah( m h “ A))<Ah * 

(13) 

= 

— ^ (^kin ^spin) • 

ZWh 

(14) 

The expansion is correct up to terms of order 1 /m^. as¬ 
suming D () i/r = O(rah), Dkif/ = 0(1) = G M . Here G ;/ 
is the gauge field, D t , the covariant derivative and we 
introduced the higher dimensional fields 

^kin(-^) 

= <AhM 1)2 'Ah(x), 

05) 

^spin(-^) 

= i// h (x ) (t ■ B(x) i^h(x), 

06) 


with 

(Tk — 2 e ii kCr ij > 7?/] * (17) 

The decoupling of the fields I s achieved by a 

Fouldy Wouthuysen-Tani (FTW) transformation (see 
El SMI) of the form 

A = exp(;X(ra)) if/ , (18) 

with 

£(w) = - 1 -D k y k + —^jkJo[Dk,D 0 \ , (19) 

2 m h 4m~ 


followed by a projection onto decoupled components 

<Ah = p+x> 'Ah = p -x> (2°) 

with 

P± = ^(1 ± 7o). P + P-= 0. (21) 

Analogous expressions for i// h and ,2^ 1 ’ are skipped 
here. Indeed, in the following we do not consider pro¬ 
cesses involving the anti-quark field, t/^, and therefore 
drop all terms containing it. 

It is worth summarizing some issues that arise in this 
formal derivation. 

• Assuming D<. = 0(1) = G ; , means that this is a 
classical derivation: in the quantum field theory 
path integral we integrate over rough fields, i.e. 
there are arbitrarily large derivatives. The renor¬ 
malization of the derived classical Fagrangian 
could then in principle result in additional terms 
of a different structure. 

• The derivation is perturbative in 1 /rah, order by or¬ 
der. This is all we want for an EFT. In this way we 
expect to obtain the asymptotic expansion in pow¬ 
ers of 1 /rah. 

• There are other ways to “derive” the form of the 
Fagrangian. One may integrate out the compo¬ 
nents P-ip , ij/P in a path integral and then per¬ 
form a formal expansion of the resulting non-local 
action for the remaining fields in terms of a series 
of local operators IfTTl . Or one may perform a hop¬ 
ping parameter expansion of the Wilson-Dirac lat¬ 
tice propagator. The leading term gives the propa¬ 
gator of the static action. 

1.5. Heavy quark symmetries 

The lowest order Fagrangian ,.S? h slat has new symme¬ 
tries. At each space-time point one may perform SU(2) 
rotations in the two-dimensional space spanned by the 
two (non-relativistic) Dirac-components. This invari¬ 
ance is the spin-symmetry, which for example exactly 
relates the correlation function of the vector current 
V^ tat = i/fy-yj. i//h to those of the time component of the 
axial current Ajl tat = 4> u ysyo |Ah- Furthermore, a phase- 
transformation i//h(x) —> e'“ (x) t/'h(T) leaves the FO ac¬ 
tion invariant. It means that the number of h-quarks is 
conserved locally at each space point. It is common to 
remove the mass term /nh^iAh from the Fagrangian. In 
Minkowski space this can be done by a time-dependent 
phase transformation of the quark fields. In the Eu¬ 
clidean it turns into an exponential factor. The physical 
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Figure 1: The lattice spacing dependence of f a( 0. 0.5), (0,0.5) and 

h( 1/4) in the quenched approximationITSi. Some data points have 
been shifted in a 2 for visibility. Different symbols refer to different 
discretisation of the static action with filled circles the original one by 
Eichten and Hill Il9l . Graph from Il8l . 


interpretation is (in both cases) that one just shifts the 
energies of the states with a single h-quark by exactly 
nib. We find it simpler to keep because a term 

of this form appears upon renormalization anyway. The 
formulation without the mass-term, however, exposes 
another symmetry, namely heavy-quark flavor symme¬ 
try which is present when more than one heavy quark 
are present. It has approximate phenomenological con¬ 
sequences. We do not need it here, mainly since in Na¬ 
ture there is no obvious partner of the b-quark which is 
heavy enough and forms bound states^] 


1.6. Theoretical status 


In section 1.3 we have emphasized the importance 
of the renormalizability of the lowest order of the ef¬ 
fective theory. In HQET we depend on the renormal¬ 
izability of the static theory. To our knowledge, this 
important property of the theory has not been proven to 
all orders in the coupling constant expansion - in con¬ 
trast to QCD. Simple power counting does not apply, 
since the static propagator does not fall off in all di¬ 
rections in momentum space. An alternative strategy is 
to prove the renormalizability after integrating out the 


4 Of course, one often considers the charm quark in an HQET ex¬ 
pansion, but on the quantitative level large corrections have to be ex¬ 
pected with a mass of around a GEV. 


static quark fields. The resulting non-local observables 
are then defined in QCD and are closely related to Wil¬ 
son loops whose renormalizability has been proven to 
all orders Il20l I2T1 l22l l23l l24ll . Indeed, in Il24ll also 
observables are considered which are very closely re¬ 
lated to the non-local observables resulting from HQET. 
Furthermore, a significant number of loop computations 
have been performed, partially to a high order ll25ft . No 
problem with the assumed renormalizability has been 
found. Also the continuum limit of the lattice theory 
has been studied in quite some detail (see e.g. lfT8l ). 
Its existence is a non-perturbative “proof” of renormal¬ 
izability. We use quotation marks, since it is a numer¬ 
ical proof only. Still, the quality of the numerical in¬ 
vestigation is very good. We demonstrate it by a graph 
from m. It shows three quantities in the static effec¬ 
tive theory, which can be computed precisely. They are 
constructed from correlation functions in a finite vol¬ 
ume with Schrodinger functional boundary conditions 
||26l 27 281 i.e. Dirichlet boundary conditions in time, 
and periodic boundary conditions in space. For the 
quark fields, the spatial boundary conditions involve a 
phase 9 via, 

i//(x + Lk) = e'V(x), (22) 

\f/{x + Lk) = e~' e t/r(x). (23) 

Such boundary conditions are (precisely speaking at 
lowest order in perturbation theory) equivalent to rais¬ 
ing the lowest momentum of the finite volume modes 
from p = 0 to p = p fl = 6 x (1,1,1) for the quarks and 
p = -pe for the anti-quarks. One can then consider var¬ 
ious types of correlation functions, see also section [4] 
Here we need ff^(9), where a relativistic quark with 
p = pa and a static antiquark with p = -p g are cre¬ 
ated at time zero and are annihilated in the middle of the 
T x L 3 space-time through the static-light axial current 
Ag tat . This correlation function is depicted on the left 
side of figure El Forming £a = f^\9)/f^ t (9') gives a 
first observable, which according to our naive applica¬ 
tion of dimensional counting and symmetries needs no 
renormalization. It is also precisely computable. For 
6 = 0,9' = 1/2 it is shown on the top of figure |T| Let¬ 
ting the quark and the static anti-quark propagate from 
the boundary at xo = 0 to xo = T, defines F\ (rightmost 
graph in figure [4] and £\. Its continuum limit is shown 
for the same kinematics in the middle graph of figure IT] 
Finally, the ratio h shown at the bottom of figure ITfis 
constructed from the boundary-to-boundary correlation 
function of a static quark-antiquark pair. For all three 
cases, the lattice spacing dependence is shown for four 
different discretizations. One can see by eye that they 
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tend to agree at a = 0. We have good numerical evi¬ 
dence for the expected renormalizability of the theory. 
Much more such evidence was seen in later works. 

In summary, we have no doubt about the renormaliz¬ 
ability of the theory. 


for a number of observables equal to the number N u of 
parameters «, present at a certain order. This procedure 
is referred to as matching. The O, do not need to be 
experimental observables. It suffices that they are ac¬ 
cessible to precise lattice simulations. 


1.7. HQET parameters 

EFTs have a number of free parameters, the w, men¬ 
tioned above. The number of parameters usually grows 
rapidly as one increases the order of the expansion. 
Concerning HQET, we provide an overview for the first 
two non-trivial orders in table [2] Depending on the 
application, i.e. the particular observable which is be¬ 
ing expanded, the number of parameters which actually 
contribute may vary significantly. For example in the 
first row of the table there are three parameters for the 
Lagrangian at NLO. These generically contribute to any 
observable, but there are exceptions. For example the 
splitting of two levels related by the spin symmetry of 
the theory (an exact symmetry at the static order) de¬ 
pends on a single parameter only, w sp j n . And for the 
computation of the spectrum of the theory, only the first 
row is needed anyway, since the other rows arise from 
the expansion of specific local fields. 

In general the growth of the numbers of parameters 
with the order is problematic for phenomenological ap¬ 
plications of EFTs. They are usually fixed from experi¬ 
ments, limiting the predictivity enormously. For HQET 
the situation is actually worse than e.g. for chiral pertur¬ 
bation theory (see j29l and references therein for a re¬ 
cent review). Even when the parameters oj, are known, 
the observables cannot be determined by perturbation 
theory in the couplings (not to be confused with the ex¬ 
pansion in l/?Hh): they are non-perturbative in the QCD 
coupling irrespective of the order in l/m^. A lattice 
QCD formulation is needed. In its absence the predic¬ 
tions of the effective theory comprise of approximate 
scaling (with powers of the quark mass) of certain re¬ 
sults from charm to bottom and there are relations be¬ 
tween matrix elements due to the enhanced symmetry 
of the static limit. Instead, with a lattice formulation, 
the theory becomes fully predictive. 

Once one has the lattice formulation, it is also natu¬ 
ral to solve the problem of the number of parameters by 
determining them from lattice-accessible quantities in¬ 
stead of from experimental ones. This is indeed possi¬ 
ble. Starting from the general idea lf30l a detailed strat¬ 
egy was developed in a number of works which we will 
discuss in the following. The basic procedure is just to 
require 

(jHQET = (pQCD J ■ = 1, . . . , , (24) 


2. Matching and mass scaling 


As simple as it is in principle, there are important 
issues to be considered and to be understood when 


eq. (24 1 is applied in practice. 


1. Observables have to be chosen such that the accu¬ 
racy of the HQET expansion is not compromised: 
energies have to be sufficiently small. There is, 
however, no point in going significantly below A ~ 
400MeV, since this scale is always present. 

2. The right hand side of eq. © has to be com¬ 
putable in lattice QCD, which means that am h <s 1 
has to be reachable. 

3. When eq. ( |24| is implemented non-perturbatively, 
it has to be imposed at a finite (the desired) quark 
mass. Since the left side is an approximation, trun¬ 
cated at a given order, the effective theory parame¬ 
ters then do depend on the matching condition im¬ 
posed. 

4. The observables <t>, have to have a good statistical 
precision in lattice computations, both in QCD and 
in HQET. 


Clearly item[2] is in conflict with section 1.2 where 
we explained that am h <s 1 cannot be reached for 
large enough volumes where finite size effects are small. 
This is avoided by having a smaller gap between in¬ 
frared and ultraviolet cutoff, defining the observables 
<I>, in a small volume. Item [I] suggests that the infrared 
momentum cutoff should be chosen around A, namely 
L = 0(1/A) = 0(|fm). 


2. 1. One-loop perturbation theory 

Before coming to a general discussion, it is instruc¬ 
tive to look at the simplest case of matching in perturba¬ 
tion theory. We consider the static effective theory. Its 
Lagrangian 


-^h tat = <Ah(' w h are + A»)^h, 


(25) 


contains a single parameter, m b 


bare 


In comparison to 


eq. (12 1 we have added a label ‘bare’ to indicate the bare 
mass parameter, fa etc are the bare fields in the regular¬ 
ized path integral. The explicit form of the heavy quark 
propagator, which we give later in eq. ( [84] ), shows that 
;;C arc drops out of all observables (at LO) except for the 
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static (LO) 


0(1 /m h ) (NLO) 

origin 

application 

Number cj. 

Number w,- 



1 

^bare 

2 

^kin ? ^spin 

j^?HQET 

all 

1 

lnZ, HQET 

2 

c Aq,i> c Aq,2 

A HQET 

A 0 

B -> tv. /b, B s -> ft. / Bs 

1 


4 


.HQET 

A k 


1 


2 


t/HQET 

v 0 

B — > tv tv. form-factors / + ,/o 

1 


4 


t/HQET 

v k 

B — > tv tv. form-factors / + ,/o 


Table 2: The number of free parameters in HQET at a given order, the specific fields where they appear (“origin”) as well as some examples for 
applications where they contribute. 


relation between the QCD quark mass and one energy 
level in the static theory, say the mass of the B-meson. 
All energy differences and all properly normalizetf]ma- 
trix elements are independent of m h f c . 

Interesting, non-trivial, matching happens for com¬ 
posite fields. We choose here the time-component of 
the axial current, 

A 0 r(x) = Z a A 0 , A 0 = if' u (x)y 5 yoi/'b(x). (26) 


For the moment the relevant property of the states 
|Q(L)) and \B(L)) is that L is the only scale apart from 
m b and A. Then, for sufficiently small L, the relevant 
QCD coupling is small and there is a perturbative ex¬ 
pansion 

M QCD (L,m b ) = M (0) + M" } cu (z)g 2 (29) 
+0(g 4 ) + 0( 1 /z ), z = Lm h . 


To distinguish it from the HQET field we label the heavy 
quark in QCD by b. The matrix element, 

<0| Aq(x = 0)|B(p - 0)) - mff b, (27) 


of the associated Hilbert space operator defines the de¬ 
cay constant / B , the only hadronic parameter determin¬ 
ing the decay rate B —> tv. For matching, it is natural to 
consider more general matrix elements 

MqcdCL. m b ) = (Q(L) | Aq(x = 0)|fi(L)>, (28) 


to define a suitable quantity <f>, (i fixed) in eq. (24 1 . In 
physical processes, L is an inverse momentum scale, but 
we will later use states in a finite periodic LxLxL torus. 
The state B(L) has the quantum numbers of a B-meson, 
while Of L) has vacuum quantum numbers. 

In generic regularizations, e.g. in dimensional or the 
Wilson lattice one, the axial current is affected by a non¬ 
trivial renormalization ,4 ;| R = When the renor¬ 

malization factor Z A is defined such that the current sat¬ 
isfies the chiral Ward identities |3ll|32 3, eq. ( [27} gives 
correctly the weak decay amplitude and thus / B . The 
current is then also scale independent. 


5 A proper mass-independent non-relativistic normalization has to 
be chosen. The standard one is (B(p')|B(p)) = 2(2zr) 3 <5(p - p'). 


in terms of renormalized coupling and mass g = g,m b . 
We will specify their renormalization scheme and scale 
when it becomes relevant. For any finite L and m b , 
the matrix elements are finite, but the large mass limit, 
m b — > oo with L fixed does not exist. It is logarithmi¬ 
cally divergent ll33l[34l . 

HcdW B ' b ~ “ (“To log(z) + Bqcd)M (0) , (30) 
To = _ l/(4zr 2 ). (31) 

This behavior has to be reproduced by the effective the¬ 
ory. As a first step, the bare static-light current 

A* tat (x) = ^ u (x)T5To^h(x) ) (32) 

is just form-identical to the relativistic one. For the clas¬ 
sical current this follows from the FTW transformation 
and beyond we just observe that there are no other di¬ 
mension three (or lower) composite fields with the same 
quantum numbers. 

Unlike the relativistic current, there are no chiral 
Ward identities which fix its renormalization. As a con¬ 
sequence the renormalized current is scale dependent. 
For example in the lattice regularization we can renor- 
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malize it by lattice minimal subtraction]^] 

< aUal (^) = ZT M (ra,go)AT\x), (34) 

Z^muat = 1 -y 0 \og(aiu)gl + O(gl). (35) 

Here a is the lattice spacing, go is the bare coupling and 
p is the renormalization scale. 

The lowest order anomalous dimension yo coincides 
with yo from the mass-scaling defined above. The ma¬ 
trix elements 

M la ‘ t (L, //) = <Q(L)|Zf l ' lat A atat |fl(L)>, (36) 

corresponding to the above QCD ones have an expan¬ 
sion in the renormalized coupling, 

AC(T,p) = M (0) +MZ(rL)g 2 + O(g 4 ), (37) 

MZirL) = (-y 0 \og(jjL) + fii at )A1 <0) . (38) 

For convenience we put the asymptotic QCD expression 
and the static one together (up to O(g 4 ,1 /z)), 

AIqcd = (1 + ( _ 7o log(ra b L) + Bqcd) g 2 )Ad (0) , 
AC = (l+(-yolog(ML) + B M )g 2 )M (0) m 

In this way one sees immediately that a finite renormal¬ 
ization of the static current 

^stat.match^; /Wb ) _ C match (m b ,ji) Ag taUat (x; ^u), (40) 

Cmatch = 1 +Ci0u/m b )g 2 + 0(g 4 ), (41) 

ci(/j/m b ) = y 0 log(p/m b ) + (Bqcd - #iat) , 


brings QCD and the static effective theory into agree¬ 
ment, 

MQCD = Cmatch AC + 0(1 /Z) (42) 

= <Q(L)| A^ tatmatch |B(L)) stat + 0(1 /Z) ■ 

We emphasise the general structure and the important 
features of the example. 


The coefficient of the logarithm in eq. (301 and the 
anomalous dimension of the current in the static 
theory, eq. ©■ match. This matching cannot be 
enforced, it is a property of the two theories. And 
it is one of the conditions for the effective theory to 
describe the asymptotics of QCD. 


6 For comparison, the renormalization factor of the relativistic lat¬ 
tice current is 


z A = i+4 1, g g + ... 


(33) 


with a pure number (no renormalization scale dependence) Z^ ! 1 , which 
can be chosen such that the chiral Ward identities hold. 


• The relative one-loop coefficient |35 36 j. 

fiQCD-fila, = -0.137(1), (43) 

is independent of the external states. They were 
chosen from two rather different classes in the two 
cited references. Again this is a necessary condi¬ 
tion for the effective theory to describe QCD. 

• Both Af (0) and Bqcd do depend on the external 
states. Let us label a matrix element for a differ¬ 
ent pair of states by just a prime. In ratios of these 
matrix elements, 

'Rqcd = A 1 q CD /A 1 qcd . ( 44 ) 

the entire renormalization and matching of the cur¬ 
rent drops out, since it is multiplicative. One then 
has effective theory predictions 


?QCD 

= Kstot + OU/z) 

(45) 

^stat 

c* Is 

+ 

s s 

ii 

g 2 + o(g 4 ), 

(46) 

.stat 


(47) 

n m 

/v stat 

(M' s tat ) (n 



(A^at) (0) 

(M S tat) (0) 



“ ^lat — B]m = 

■ Bqca , _ Bqcd ■ 

(48) 


• The particular number of eq. ( |43| depends on the 
renormalization scheme for the static current. Here 
we chose minimal subtraction, a scheme which is 
not independent of the regularization. Therefore, 
Bqcd - B\ m depends on the details of the regular¬ 
ization chosen in lf35l [36 l. It is valid for the O (a) 
improved Wilson lattice regularization. 

• Of course the matrix elements of the matched static 
current in eq. ( |42| do not depend on any details 
of the regularization. Their finite renormalization 
has been chosen to match QCD. This is unique. 
In eq. (|42| we did not indicate the one-loop na¬ 
ture. Indeed, we expect this equation to hold to 
all orders in the coupling and also beyond, non- 
perturbatively. 

• We have nowhere given the renormalization 
scale/scheme for coupling and mass. At the one- 
loop order all expressions are independent of it. 
The scheme only matters for the renormalization 
factor of the current itself. In the following section 
we will also discuss a convenient choice of renor¬ 
malization scales. 
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2.2. Higher orders in the coupling and renormalization 

group invariants 

In this section we explain on the one hand what is 
known at higher orders in the coupling and on the other 
hand, how one passes to renormalization group invari¬ 
ants, which are independent of schemes and scales. In 
particular they allow for a clean factorization of observ¬ 
ables into a non-perturbative matrix element and a mul¬ 
tiplicative matching function, which has a perturbative 
expansion. This separation makes efficient use of the 
high order perturbative information accumulated over 
the years ll?3l l34l l37l l38l l39l l40l l4ll l42l l43l l44l l45l. 
We note, however, that we only know how to apply this 
strategy to the lowest order, static, effective theory. A 
hurried reader may therefore skip this section and pro¬ 
ceed to the following one. 

2.2.1. RG functions and invariants 

It is well known, that fixed order perturbation theory, 
where all observables are expressed in terms of a renor¬ 
malized coupling and masses at one fixed renormaliza¬ 
tion scale is not the best choice. In fact, if scales rather 
different from the renormalization scale are relevant in 
the observables, large logarithms multiplying powers of 
the coupling are present and the accuracy of perturba¬ 
tion theory is not good. This is a reason to consider 
running coupling and mass, i.e. coupling and mass as 
functions of the renormalization scale, 

g(p), Thi(p ), (49) 

and the associated renormalization group equations 
(RGE). Here i runs over the different flavors and all 
masses are defined in QCD, also the mass of the quark 
treated by HQET. A consequent application of the RGEs 
is to use their solutions and express all observables in 
terms of renormalization group invariants. 

Here we do not give an introduction to the RG but 
just recommend Ref. l46l . We mainly describe the rel¬ 
evant formulae in order to apply them to our matching 


problem. We work in an unspecified massless renormal¬ 
ization scheme, where the renormalization factors do 
not depend on the masses. Consequently the renormal¬ 
ization group functions do not depend on the masses. 
Examples of massless schemes are (modified) minimal 
subtraction in dimensional regularization (MS, MS) or 
lattice regularization (lat) or a Schrodinger functional 
scheme (SF) 123147). The latter is independent of the 
regularization. 

Our renormalization group (RG) functions are de¬ 
fined through 


= P(g) 


H dnij 
nil dpi 


T(g) , 


r dM l,* 

M lm & 


= y(g) ■ 


(50) 

(51) 


Apart from the running coupling and running quark 
mass we here consider the matrix element Af s ' tat of a 
(multiplicatively renormalizable) composite field renor¬ 
malized at scale p. We label it with a superscript s for 
the scheme to remain consistent with previous notation. 
One may identify it with eq. (36 1 , but other matrix el¬ 
ements of general composite operators are possible as 
well. The RG functions have asymptotic expansions 


m 

S ~° ~g 3 {^o + g 2 bi + ...} , 

(52) 


1 38 

bi = , (102 N { ) , 

1 (47t) 4 ' 3 ' 

(53) 

t (£) 

g ~° -g 2 {<7o+g 2 ^i +, 

(54) 


d 0 = 8/(4tt) 2 , 

(55) 

7(g) 

g ~° -g 2 {ro + g 2 y i + ...}. 

(56) 


The integration constants of the solutions to the RGE 
define the RG invariants 


A = ( fig(g)l l = Li(bog 2 ) /l/Uo, e 1/(2i °s 2) exp j- 
M, =<p m (g)mi = Thj (2bog 2 )~ d ° /2b ° exp j - dv 
ACt 1 = Vstat® Ktat = Ml tat \2b o g 2 Y y "' 2b0 


dv 


1 1 _ bi 

P(x) + b 0 x 3 b\x 


t(,t) _ do_ 

P(x) b 0 g 

y(x) _ jo_ 
P(x) b 0 x 


(57) 


-X 


dv 


r 
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where g = g(p), m t = m, ,(p), M* tat s Al s s tat (p). There 
are no corrections to eqs. ( |57j ) with exact RG functions. 
Instead, when we only have access to their perturbative 
expansions such as eq. ( [52] ) there are truncation errors 
in the functions <f,(g). For example, in order to ob¬ 
tain a numerical result for A given a pair of numbers 


and we have 
AIqcd 


Cmatch(g*) X Ad s t a t(/r) 
_ Cmatch(g*) y^jRGI 


0(1 /m h ) 

stat + 0(l/w h ). (66) 


for g(p), p, the integrand [... ] in eq. (57 1 may in princi¬ 
ple be expanded as [... ] = C 2 X + c^x^ + 0(x 5 ). Knowing 
C 2 , the error made is (asymptotically) ~ C3g 4 /4 and the 
same is true when we insert the truncation of / 3(x ) and 
then integrate. When we discuss numbers, we will al¬ 
ways expand the RG functions, not the full integrands. 

Physical observables can be written as 


(^stat(g*) 

The g*-dependence is equivalently to the mass depen 
dence and defines another RG function y matc h, 

* SA-Iqcd 


yrnatch(W) — 
with asymptotics 


AIqcd dm. 


(67) 


V = P(KAM l \,{p l \) 


(58) 


making it manifest that they do not depend on a renor¬ 
malization scale p, 


H-r-P = 0. 

dp 


(59) 


In the following section we express the lowest order 


yma,ch(g) ~ - yog - rr V + ■ • • ■ (68) 

Because the leading order coefficient yo is unchanged, 
“match” is just another renormalization scheme for the 
axial current Ag tat . In this scheme, at scale p = ;«», its 
matrix elements are equal to the QCD ones up to order 
1/mh. We therefore refer to it as the matching scheme. 

Finally, a complete transition to renormalization 
group invariants is achieved by changing to the RGI 
mass scaling function 

M <9C PS 


HQET approximation of matrix elements AIqcd in this 

rWM/A) - M 

way and define a RGI mass scaling function, which 

Prs(M/A) _ m 


C PS dM 


(69) 


describes the mass-dependence beyond the asymptotic 
form eq. 0 

2.2.2. Mass scaling 

As a first step we choose a suitable renormalization 
scale, 

p = m, , (60) 

where the solution of 

m,-m(m,), (61) 

defines m,. We further use the shorthand 

g, = g(m,) . (62) 

The coupling g, can be determined for any value of 
M/A by combining the first two equations eq. © to 

A _ <P g (g») 

M g} m (g„) 


ymatch (g*) 
1 - T(g„) 


and to 


Mqcd = Cps(M/A)xMf° I + 0(l/m h ) 
C PS (M/A) = Cpnatch^*) / ^stat(,?*) 

g* , Tmatch^) 1 


exp 


/ 


ix ^wr)‘ ,70) 


Everywhere it is understood that g, = g,(M/ A). In the 
last equation, the integral at the lower bound is to be 
understood as in eq. ©• 

At leading order in 1 //n h the conversion function Cps 
contains the full (logarithmic) mass-dependence, while 
the non-perturbative effective theory matrix elements, 
are mass independent numbers. 

An interesting application is the asymptotics of the 
decay constant of a heavy-light pseudo-scalar (e.g. B)0 


= exp 


-f 


■*•("/ A) 1 - t(x) 
ax - 

m 


r co [ln(M/A)] r °^ 2i ° aRGI 

rps ~ -;--At, 


(63) 


^Jmps 

X[1 + 0([ln(M/A)] _1 )]. 


(71) 


The solution of this equation defines a function 
g*(M/A). 

With the above choices the matching function simpli¬ 
fies to 

Gnatch dll. - ft!. ) = G ma t c h(g*) (64) 

= 1+ Cl (i)g; + ..., (65) 


In perturbation theory, the logarithmic corrections are 
computed by solving eq. ( |63] for g, and then integrating 
eq. ( [70] ). 

7 Note the slow, logarithmic, decrease of the corrections in eq. ED- 
We will see below, in the discussion of figure|2] that the perturbative 
evaluation of Cps(Mb/A) is somewhat problematic. 
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1-loop - 2-loop 3-loop | | . 2-loop - 3-loop 4-loop | 




Figure 2 : The function y m atch(g) as a function of g 2 = g 2 for Nf = 3 flavors in the MS scheme. On the left we show y ma tch = y^° tch for the time 
component of the axial current, on the right we show the difference y^ atch - 7^° ltch - Note that at 1-loop order the latter vanishes. 


2.3. On the accuracy of perturbation theory 

When one evaluates functions such as Cps in a given 
order of perturbation theory, various quantities enter 
such as the beta-function, the quark mass anomalous di¬ 
mension. Apart from y matC h, these all have a well be¬ 
haved perturbative expansion in the MS scheme, see ap¬ 
pendix A.2.2 of 0 for a table of the coefficients. Keep¬ 
ing this in mind, we just discuss y matc h- In the left graph 
in figure [2] we plot different orders of y nlatc h ■ For the 
larger values gl in the plot one may get worried about 
neglecting higher order terms. Note that gl is around 
2.5 for the b-quark and it is out of the range of the graph 
for the charm quark. 

However, a more serious reason for concern derives 
from the right hand side graph. There the difference 
of the anomalous dimensions for 14 and ,4 (l is shown. 
For such differences perturbation theory is known to 
one loop higher Il48l and the perturbative coefficients do 
grow further. Asymptotic convergence seems to be use¬ 
ful only for rather small couplings or masses far above 
the b-quark mass. At the b-quark mass every known 
perturbative order contributes about an equal amount. 
Since we do not understand the reason for this behavior, 
it raises concern about using perturbation theory for the 
matching functions. 

Let us emphasize, that the bad behavior is easily 
traced back to the function C matc h and was noted in 148 1. 
We tried earlier 0 to rearrange the perturbative series 
in order to find a more stable perturbative prediction, but 
we did not succeed. 


2.4. Matching at NLO in \ / 

For quantitative phenomenological results one has to 
compute also 1 /mh corrections in HQET. Is it consistent 
to match perturbatively as we discussed in the previous 
sections? We saw that the uncertainty due to a trunca¬ 
tion of the perturbative matching expressions at /-loop 
order corresponds to a relative error 


A(Cps) 

Cps 


A(A 4 qcd) 

WIqcd 


L? 2 ('m«)]' 


2£>o ln(m*/AQCD) 


(72) 


As m t is made large, this perturbative error decreases 
only logarithmically. It becomes dominant over the 
power correction which one wants to include by pushing 
the HQET expansion to NLO, 


A(Cps) m,3> A A 

-7 - » — 

Lps ill * 


(73) 


With a perturbative matching function, one does not per¬ 
form a consistent NLO expansion such that errors de¬ 
crease as 1/m”. 

A practically even more serious issue is that at NLO 
one has to deal with the mixing of operators with lower 
dimensional ones. For example 0^ n = (// h l.)A// h mixes 
with i/r h Z)()i^h and i/'h'Ah- In this situation mixing co¬ 
efficients are power divergent ~ a ". In the example 
we have n = 1,2. Subtracting power divergences in 
perturbation theory and then computing the matrix ele¬ 
ments non-perturbatively always leaves a divergent re- 
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mainder. The non-perturbative continuum limit of ma¬ 
trix elements of perturbatively subtracted operators does 
not exist. 

We are lead to conclude that it is necessary to per¬ 
form matching and renormalization non-perturbatively. 
The only alternative is to supplement the theory by as¬ 
sumptions. Namely one may assume that at the lattice 
spacings available in practice, power divergences of the 
form gpV (m b a) n are small since m b a > 1. This then has 
to be combined with the assumption that the b-quark 
is not large enough to be in the asymptotic region of 
eq. {73]). 


2.5. Non-perturbative matching 


2.5.1. Scope 

A true non-perturbative matching as in eq. (24 1 elim¬ 
inates the problems we just discussed: slow asymptotic 
convergence of PT, power divergent remainder terms in 
the subtraction of lower dimensional operators. We do, 
however, still have to cope with the condition m b a <sc 1 
for the computation of (D qcd . Before we explain the 
strategy to deal with this, we would like to emphasize 
a point about the magnitude of 1 /wh corrections which 
sometimes leads to confusions. 

The size of these corrections does depend on the 
matching conditions used: after imposing a generic con¬ 
dition, eq. ( |24| ), the decay constant of the B-meson once 
computed at LO and once computed at NLO in l/m b 
differ by an amount Of 1 /m b ). However, how much this 
difference is exactly, depends on which set {<!>,} was 
chosen for matching the theories. We are free to choose 
<f>i = me , d >2 = (0| Ay|B(p = 0)). With this choice, 

the decay constant of the B-meson is exact at all orders 
in 1/wh. 


This rather trivial fact has to be remembered. The ex¬ 
act size of the corrections is only defined, once one has 
fixed how the matching is performed. What we said be¬ 
fore implies that one has to define non-perturbatively 
how the matching is performed. Only then does the 
splitting of a HQET result into different orders acquire 
a precise meaning. 

Let us illustrate the consequences on a frequently dis¬ 
cussed example, namely the mass formulae 


m 


av 

B 


Awb 


i[m B +3m B *] (74) 

m b + A + —/li + 0(1 /ml) (75) 

2m b 

2 , 

m B > - m B =- A 2 + 0(1 /ml) (76) 

m b 


with (ignoring renormalization) 


M = < B\OUB >, A 2 = ^(BIO spm IB). (77) 

The quantity A is referred to as “static binding energy” 
and Ti as the kinetic energy of the b-quark inside the 
B-meson. Also here, depending on how one formulates 
the matching conditions, one changes A by a term of 
order Aqcd- Similarly, the kinetic term 4i/(2;nb) has a 
non-perturbative matching scheme dependence of order 
A 2 QCD /m b and thus A t itself has a matching scheme de¬ 
pendence of order A ( 2 )C| ) . In fact we could set Oi = m™. 
At static order this means m b + A = ;;f' v and at NLO it 


is eq. (|75|) without the 0(1 /mi) correction. If both are 


to be valid, one has Ai = 0. 

Lor this reason, our scope is not to compute (and 
therefore first define) quantities such as A but to com¬ 
pute physical observables such that they are correct up 
to corrections of order 1 //;?-. Lor this purpose we need 
to determine the bare parameters to,-. We do not need to 
define renormalized ones. We ask the reader to keep in 
mind that since a>, are bare parameters, they depend on 
the bare gauge coupling (equivalently the lattice spac¬ 
ing) and the heavy quark mass, where we may choose 
the RGI mass, M. Lurther dependences on the details of 
the discretisation are kept implicit. 


2.5.2. Strategy 

The general idea how to satisfy m b a <c 1 is to per¬ 
form the matching step in a finite space-time volume of 
linear dimension L = L\ |[30j| , exactly in line with the 
idea how to cover the scale hierarchy in the computa¬ 
tion of running couplings l50l . Choosing L/ 1 ~ A, the 
HQET expansion of correlation functions at distances 
|x| = O(Li) is valid. L\ ~ 0.5 fm is a reasonable choice. 
This allows for lattice spacings of a - 0.02 fm and sig¬ 
nificantly smaller. One can then perform continuum ex¬ 
trapolations with several points satisfying m b a <1/2 
and determine 

O qcd (Li , M, 0) = lim ® qcd (L 1 ,M,a), (78) 

a-> 0 

as indicated on the left side of figure|3] One is now free 
to use the matching condition, eq. ( |24| t, to determine the 
COi at a given resolution L\/a. Explicitly, we write the 
HQET expansion in a matrix notation 

4> HQET = rt + cpu, (79) 

with a Nu x N w matrix cp and a homogeneous static part 
rp, i = 1,..., N 0J . Both rj and ip can be computed by nu¬ 
merical simulations of HQET. Solving eq. {79| defines 
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Figure 3: The ALPHA collaboration strategy for matching of QCD and HQET. The bottom row denotes the continuum limit obtained by an 
extrapolation of the results in each column. Graph from ED. 


a first set of HQET parameters 

u(M, a ) = (80) 

V~ l iU, a)( cD qcd (L, , M, 0) - r/(L\, a)) . 

For reasonable resolutions, say L/a > 8, where dis¬ 
cretisation errors may be assumed to be small, we have 
a < Li/S w 0.06fm. Remembering that one should al¬ 
ways have a range of lattice spacings in order to carry 
out a continuum extrapolation, we realize that such lat¬ 
tice spacings are still too small to perform the large vol¬ 
ume HQET computations to determine the physical en¬ 
ergy levels and matrix elements. 

For this reason, the full strategy includes a so called 
step scaling step, to reach the same cD HQET but for L = 
L 2 = 2L\\ 

® HQET (L 2 ,M,0) (81) 

= lim r/(L 2 , a ) + ip(L 2 , a) u>(M, a) . 
a—>0 L J 

This set of observables then serves to determine the de¬ 
sired parameters 

a>(M, a) = (82) 

y _1 (L 2 , a)(<t> HQET (L 2 , M, 0) - i](L 2 ,a)) . 


Here lattice spacings can be used which are suitable for 
large volume HQET computations. The overall strategy 
is depicted in figure [3] It is explained in more detail in 

m. 


2.5.5. The electroweak heavy-light currents 

For the phenomenological treatment of (semi-) lep- 
tonic B-decays in the standard model, one needs both 
the axial current and the vector current. Our discus¬ 
sion above has to be generalized accordingly. The fi¬ 
nite renormalization (for matching) of space and time- 
components is different, but in a minimal subtraction 
scheme, all currents have the same anomalous dimen¬ 
sion in the static effective theory because they are re¬ 
lated by spin symmetry and the chiral symmetry which 
emerges when the light quarks are massless. Details can 
be found in Q. 

Non-perturbative matching conditions for the full 
set of currents have been discussed in GQ. They 
have partially been investigated in perturbation theory 
11521 37 54l|. see 1331 for more details. 
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3. Discretized HQET 


We need a lattice formulation in order to treat HQET 
beyond perturbation theory in the QCD coupling. The 
static action discretized on a hyper-cubic lattice is 

-2p - t -^ <AhW[Vo + mjTMhM, (83) 

with Vp the gauge covariant backward derivative. Com¬ 
pared to the form written down first by Eichten and Hill 
HI we just added the mass term. The static propagator 
in a gauge background is 


G h (x,y) = 


with 


8(xq - yo) S(x - y) e 
x P(y, x; 0) f P + 


o-yo) 


„bare _ 


- ln(l + am bdre ) , 
a 


(84) 

(85) 


where P(x, y; 0) parallel transports fields in the funda¬ 
mental representation from y to x along a time-like path 
and the lattice 0-function is 0(f) = 0 for t < 0 and 
0(f) = 1 for t > 0. It is a simple exercise to obtain 


eq. (84 1 from the defining equation 


(V5 + < re )G h (x,y) = <5(x-y), 


1 + am bme 
h 


where Vq is the covariant derivative with respect to the 


argument x. The propagator eq. (84 1 is valid in any 


gauge field background and therefore shows that all 
correlation functions with a static quark have a mass- 
dependence C(x) = C(x)\ lr r=0 xe-< e '*°l. Equivalently 
we see that the mass m^ dre in the Lagrangian just yields 
an energy shift 


E qcd = E s ' 


n^=0 + m h 


bare 


( 86 ) 


We now briefly discuss the most important features 
of the discretized static effective theory. The action is 
(non-perturbatively) O (a) improved without adding any 
counter-terms and therefore without tuning any param¬ 
eters. This is a consequence of the local conservation 
of the b-quark number and the spin symmetry, which 


are exactly preserved in the formulation eq. (83 1 on the 
lattice J56j. These symmetries do not allow for di¬ 
mension five fields (containing 0"h, <Ah) ' n Symanzik’s 
effective Lagrangian describing the discretisation er¬ 
rors. As a consequence energies scale to the continuum 
with lattice spacing errors proportional to a 2 log(aA) ,; at 
small lattice spacing Improvement terms of the form 


a60, [ 60\ = 4 are, however, necessary for fields O such 
as O — A? at . For the matrix element determining /b, 
there is a single such counter-term l56l . After includ¬ 
ing it with a properly chosen coefficient, one again has 
()(a 2 ) scaling to the continuum limit. 

The above features are very attractive. There is, how¬ 
ever, a (numerical, Monte Carlo) issue with the static 
theory, eq. ( [83] >. It has its origin in the fact that 

£ Stat |„,bare =0 ~ ( - ^ + 0(fl°)) + O(^) , (87) 

has a linear divergence (which is then cancelled by 
m bdrc ~ -g 2 / l ^/a in the physical energies). The nu¬ 
merical value is 

r (1) = 0.1685 for the action eq. ( |83j ). (88) 

The problem is not - as originally was thought - to de¬ 
termine m bdle non-perturbatively. Rather the statistical 
errors of, for example, two-point correlation fuctions of 
a static B-meson with m bdrc = 0 behave approximately 
as^EUEO) 

stat. error x °~“^l N e“ m ''- r,l/2 (89) 

while (at least at small a) the correlation function scales 
as 


correlator ~ . (90) 

We then have 

S1 S nal _ 54.S c -{E"*-mJ2]x ri (- 91 ) 

noise 661 s 

where as before E sidt for m^ se = 0 enters. While 
methods to use translation invariance with the help of 
“stochastic sources” ED can reduce the pre-factor 66 In, 
the correlators inevitably disappear in the noise at some 
time xq = 0(a/r (1) ). We here set go ~ 1 as appropriate 
for the considered gauge action. Inserting further that 
small lattice spacings means a < 0.05fm, we see that 
it is very difficult to have reasonable precision around 
and beyond one fm distance. However, this is the time 
separation where one can be confident that the desired 
ground state dominates. This problem brought progress 
in the static theory (and therefore in HQET) to a stop in 
the beginning of the nineties. 

Our description of the issue already suggests the 
solution. The coefficient r (<> is regularization depen¬ 
dent. Acceptable alternative discretizations can be 


8 See [57:;(581 for a discussion of the logarithmic modifications. 


^Remember that we can set m!’ are = 0 due to eq. I86>. 
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found which have much smaller values mm, for ex¬ 
ample 

r (1) = 0.0352 for “HYP2” action HU . (92) 

Our previous estimate is now pushed to time separations 
of around 4 fm. In practice it turns out that the estimate 
is overly optimistic, since it neglects the finite terms in 
E stat . The effect of the size of r ( V) is also seen in figurefy 
where the results with the Eichten Hill action, eq. (pf3p7 
have much larger statistical errors and these do grow 
towards the continuum limit, where the time separation 
in the correlation functions diverges when it is measured 
in lattice units. 

It is interesting to compare the HQET situation 
(HYP2) to the one for Nucleon matrix elements. The 
energy difference f? stat - m„/2 has to be compared to 
'«nucieon - 3m^/2. At a pion mass of m„ = 350MeV 
and a lattice with a small lattice spacing, a = 0.05 fm 
we compare about 900 MeV for the former to 600MeV 
for the nucleon. For HQET this gets (slowly) worse 
for even smaller lattice spacings and for the nucleon 
it gets quickly worse towards smaller pion masses. In 
both cases it is an unpleasant mass scale that governs 
the time-dependence of the signal-to-noise ratio. 

Given that correlation functions at 2-3 fm are not re¬ 
ally accessible, one would like to have methods which 
work reliably at smaller distances. A technique that ac¬ 
celerates the approach of masses and matrix elements 
to the ground state quantities is the use of correlation 
function matrices constructed from a few smeared in¬ 
terpolating fields |63. 64|. These correlation functions 
can be analysed by the “GEVP” method lf65l, 1661 l67l 
and then a systematic acceleration of the (asymptotic) 
convergence to the ground state is present. Physically it 
means that one uses trial wave functions in a variational 
calculation. The method is easily modified to be appli¬ 
cable also to matrix elements 1671 [681 and to effective 
theories such as HQET. We note that the acceleration 
of the asymptotic convergence is proven mathematically 
E2, but it is less clear where asymptotia sets in in prac¬ 
tice. The numerical results of section0are based on the 
GEVP technique. 


4. Verification of HQET 

In a first step in non-perturbative investigations of 
HQET, it is of considerable interest to verify the mass¬ 
scaling of QCD observables predicted by the lowest or¬ 
der effective theory. This can only be done with the 
help of lattice gauge theory, where we can change the 
quark masses. As for the matching strategy discussed 


above, a natural choice is to focus on observables in 
a finite volume. For reasons to be explained shortly, 
Schrodinger functional boundary conditions were cho¬ 
sen with a TxLxLxL geometry, where periodic bound¬ 
ary conditions (up to a phase 8 for the quarks, see sec¬ 
tion 1. 61 are present in the LxLxL space. Without dy¬ 


namical fermions the verification was carried out with 
L » 0.2 fm T - L 1691 and with N{ = 2 dynamical 
fermions with L « 0.5 fm, T — L f70| . 

We here focus on Nf — 2 observables constructed 
from the correlation functions 


/ao(*o>0) = 

~ J d 3 yd 3 z (A 0 (x)4(y)y 5 £(z)) , (93) 

kj}(xo,0) = 

4 2 / dWz (\4(x)4(y)y^(z)) , (94) 

as well as the boundary-to-boundary correlations 

F ' = -34/ d5 ^ d VdV 

(iXfyYs&V) 4(y)75<Ti(z)) , (95) 

* = - 44 / dJyd!,dVdV 

(^(y')r^b(z') 4(y)rt4i( z )) • (96) 

They are illustrated in figure [d] Fields 4 X can be 
thought of as quark fields at the boundary xq = 0 while 
4,4 are located at xo = T. The above correlation func¬ 
tions are gauge invariant despite the different locations 
of 4 and 4 This is due to the fixed gauge fields at the 
boundary; it is a very useful feature of the Schrodinger 
functional which allows for the projection onto quark 
momenta 


p = p e = (0,e,0)/L (97) 

in the above correlators (and p = -pg for the anti¬ 
quarks). 

As an example we discuss the function /a o (xo,60 
in some detail. It describes the creation of a (finite- 
volume) p = 0 heavy-light pseudoscalar meson state, 
|y> B (L)), through quark and antiquark boundary fields 
which are separately projected onto momenta p = ±pg. 
This state “propagates” an interval xo in Euclidean time. 
From the upper boundary a state with vacuum quantum 
numbers propagates a time-distance T - xo. The corre¬ 
lation function /a 0 (T /2, 8) can thus be written in terms 



Rainer Sommer 


(b ^5 / 


16 


T 


^o(*o) 



V k (x o) 





4 rs 


Figure 4: The Schrodinger functional correlation functions /a, and /i. For K \, in the rightmost diagram 75 is replaced by y*. 


of Hilbert space matrix elements. 


/ Ao (772,0) = Z-^fiCDIAoIBCD), (98) 

TO) = e-™ /2 |^ B a)>, (99) 

|Q(L)> = e- ra/2 |^ 0 a)). (100) 

where H is the QCD Hamiltonian and 

Z = <0(L)|Q(L)). (101) 


W.L)) denotes the Schrodinger functional intrinsic 
boundary state. It has the quantum numbers of the vac¬ 
uum. All states appearing in our analysis are eigenstates 
of total spatial momentum with eigenvalue zero. Their 
dependence on 6 has been suppressed. The time evolu¬ 
tion operator e suppresses high-energy states ex¬ 
ponentially. When expanded in terms of eigenstates of 
the Hamiltonian, |Q(L)) and | B(L)) are thus dominated 
by contributions with energies of at most A E = 0( 1 /L) 
above the ground state energy in the respective channel 
(recall that we take T = O (L)). 

This explains why, at large time separation xq » 
1 /;«h, HQET is expected to describe the large-mass be¬ 
havior of the correlation function, also in the somewhat 
unfamiliar framework of the Schrodinger functional p*| 



1/z 


Figure 5: Extrapolations of ips/ps(#i, 62 ) and Tv/v(#i> 02 ) to *h e 
M —> 00 limit for three combinations of { 61 , 62 ). The dimension¬ 
less heavy quark mass z = ML is used as a variable. Black circles 
show the continuum results of the corresponding quantity computed 
in the static approximation. All data points were first extrapolated to 
the continuum limit ESI- Three different combinations of 6 \, 62 are 
shown, where 0 < 0/ < 1. Graph provided by P. Fritzsch, based on 

Go). 


Equations similar to the above hold for one only 
needs to replace pseudoscalar states by vector ones. 
Finally, the boundary-to-boundary correlator is repre¬ 
sented as 


F\ = Z~ l (B(L)\B(L)). (102) 

Since the boundary quark fields £,£,... are multiplica- 
tively renormalizable |28), this holds also for the states 
\ipo(L)) and |^b(L)>. 


10 More generally, HQET applies to correlation functions at large 
Euclidean separations. 


It now follows that the ratios 


Eps(0) 

„ / Ao (r/2,0) 

V4(0) 

(103) 

Yv(ff) 

k^T/ 2 , 6 ) 

— - - , 

(104) 

Eps/ps(#i, O 2 ) 

Tps( 0 i) 

r PS (0 2 ) ’ 

(105) 

Y\/\(d 1, 62 ) 

E v (0i) 

Yv(9 2) 

(106) 

are finite quantities. 

As is immediately clear from the 
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foregoing discussion. 


Y PS (L,M) 


<Q(L)| Aq| B(L)) 

|| |n(L)> || || | B(L)) || 


(107) 


(or Y\:(L, M )) becomes proportional to the pseu¬ 
doscalar (or vector) heavy-light decay constant as L —> 
oo. At any fixed L, the large-M behavior of these quan¬ 
tities is indeed described by HQET if it is the correct 
effective theory. 

Numerical tests of this equivalence were originally 
done in the quenched approximation j69l . with L ~ 
0.2 fm, profiting from the full improvement of the the¬ 
ory with mass-dependent improvement terms im Here 
we show more recent results with Nf = 2 dynamical 
fermions. Based on the simulations |[49l the following 
results are due to Fritzsch, Garron and Heitger l70l . 

The first two HQET predictions are 


Jps/ps(0t,02) = E s tat(0i,02) + O(l/M) (108) 

YvM0 u e 2 ) = T stat (0 1; 0 2 ) + O(l/M). (109) 


Due to spin symmetry, the static limit E s t a t(fli, 0 2 ) is the 
same for the vector channel and the pseudoscalar chan¬ 
nel. It can be computed by replacing the relativistic 
fields by the static ones. The comparison is shown in 
figure [5] extrapolating the relativistic results with just 
a linear function in 1/z for z = ML > 5. Results of 
the extrapolation are in very good agreement with the 
numbers obtained directly in the static approximation, 
demonstrating at the same time the correctness of the 
effective theory and the usefulness of such observables 
for matching HQET to QCD. 

Two more predictions read 


ips(fl) 

Cps(M/A) 

Yy(0) 

C V (M/A) 


X^(ff) + 0(1/M), 
X*®(0) + O (1/M). 


(HO) 


The matrix element X^® is defined as Tps, but with the 
static quark field and it is the renormalization group in¬ 
variant one, see eq. (57 1 . For the time-component of 
the axial current, the factor i^statC!) in that equation and 
the renormalisation factor Z)!‘ ll " SI were determined non- 
perturbatively in a Schrodinger functional renormaliza¬ 
tion scheme l72l . Thus X^® is known without pertur¬ 
bative uncertainties. 

The comparison is shown in figure [6] The agreement 
between the extrapolation of relativistic results and the 
static effective theory is not as convincing as in figure[5] 
In fact, for Tps there is a disagreement, but only at the 
level of 1-2 sigma. This may therefore be a statistical 



Figure 6: Comparison of static extrapolations of Yps(0) and Yy(0) 
to the non-perturbative HQET results in the continuum (filled circle). 
Conversion functions Cps,Cv evaluated at 3-loop order of perturba¬ 
tion theory are used. The extrapolation, quadratic in 1/z, uses data 
with 1 /z < 0.2 only. Graph provided by P. Fritzsch, based on Go). 


effect, but it may also be due to the perturbative approxi¬ 
mation of Cps, which is somewhat doubtful as discussed 
earlier. In comparison, in figure [5] the perturbative fac¬ 
tors Cps, Cy drop out. 

We have only shown two out of a number of tests. 
Considering them all [69, 70|, we conclude that the ef¬ 
fective theory is very well tested, but it is safer not to use 
conversion functions Cps, Cv from perturbation theory, 
even if they are determined at three-loops, i.e. with rel¬ 
ative errors of order g 6 (m *). 


5. Numerical Simulations and results 

We now turn to a discussion of numerical results skip¬ 
ping most of the details of the Monte Carlo simulations. 
These separate into two categories. One part concerns 
small volume simulations L < lfm with Schrodinger 
functional boundary conditions. These have been dis¬ 
cussed in [73] ;74j [73, 76 77j. Many such simulations 
had to be carried out in order to control the renormal¬ 
ization and matching. They required tuning of the bare 
parameters such that continuum limits can be taken at a 
fixed volume in physical units and at vanishing dynam¬ 
ical quark mass. 

A second part is then necessary in large volume, with 
L > 2fm, nijrL > 4. These simulations are rather uni¬ 
versally useful and also much more expensive in the nu¬ 
merical effort. Hence they have been carried out in co¬ 
ordination with several European groups, by the CLS 
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effort 178], see in particular Da. The simulations were 
carried out down to pion masses of 200MeV on lattices 
with up to 128 x 64 3 points. They were possible due to 
a significant improvement of algorithms and their im¬ 
plementations lt80l 18T11821 [83] [84l [85. [86j| starting from 
the principle of the HMC algorithm 1871 l88ll . For re¬ 
cent reviews covering lattice QCD algorithms we refer 
to 1 8911901 . Some more details can be found in l9ll . 

All HQET computations and strategies were devel¬ 
oped and tested in the quenched approximation, i.e. 
QCD with all valence quarks but a vanishing number 
of sea quark flavors: Nf = 0. We will mention these re¬ 
sults only for comparison, while we discuss the results 
with Nf = 2 quark flavors in more detail. A dynamical 
strange quark has not yet been included in the HQET 
simulations of the ALPHA collaboration. 

5.1. The B*B n coupling. 

As explained above, the available large volume simu¬ 
lations still have unphysical quark masses and results 
need to be extrapolated to the point where the light 
quark masses are the physical ones. The natural way 
to carry out such an extrapolation is with the help of a 
systematic expansion in the light quark mass. Again this 
means that one uses predictions from an effective field 
theory which implements the expansion. It is heavy me¬ 
son chiral perturbation theory (HMChPT). The fields in 
this effective theory are a triplet of pion fields as well 
as a (static) B-meson and a B* meson field. The expan¬ 
sion is a combined expansion in l/m^, in the squared 
pion momenta and in the light quark mass, each count¬ 
ing as one power of the expansion variable, Aqcd//uh = 
0(p 2 /A^ CD ) = 0(ot up /Aqcd). where as before we work 
in the B-meson rest framef^l At lowest order, the ef¬ 
fective theory contains five parameters, usually called 
low energy constants. They all refer to the chiral limit 
and ;;i h —> oo. There are the pion decay constant, /, 
the light-quark condensate, the mass of the B-meson 
and the three-point coupling coupling B* B n denoted 
by g mm [981- Apart from g the couplings can be 
taken from experiment (the light quark condensate is 
removed from the list by taking mi instead of m up as 
mass-parameter). Typical predictions of HMChPT are 


(with £ = w 2 /(8tt 2 / 2 )) 

m B = mf r + f„ (3 V2t ?g 2 ? 12 + + 0(£ 2 )) 

( 111 ) 


and 




( 112 ) 


3 1 + 3g 2 

4 2 


£l°g(£) + cr f ^ + 0(£ 2 ) 


The coefficients of t he lea ding non-a nalyt ic terms, 
£ 3 / 2 and ^ log(£), in eq. (111 1 and in eq. (1121, respec¬ 


tively, are given in terms of g. This coupling is there¬ 
fore better determined before using the HMChPT for¬ 
mulae to extrapolate from unphysical quark masses to 
the physical one. 

A determination of g just means that HMChPT is 
matched to HQET at the lowest order in l/mt,; at this 
level HQET is regarded as the fundamental theory. The 
matching condition can be written to directly give g via 


8 = -<B u (0)|Af(0)|B*. + (0)>, 


A p u (x) = 


(113) 

(114) 


where i^d(iAu) annihilates a down(up) quark and the in¬ 
dex k — 1, 2, 3 is not summed over. The non-relativistic 
normalization of states given earlier is used; in finite 
volume it is (B () (p)|B 0 (p)> = <B*(p)|B*(p)> = 2L 3 = 2V, 
where L is the linear size of the torus. 

Profiting from a newly developed method for its com¬ 
putation f68l . the matrix element eq. ( | 1 14| >, was deter¬ 
mined in j92l with a much better precision than it was 
possible before. A comparison to other computations is 
shown in figure [7]for various pion masses (i.e. dynami¬ 
cal quark masses) and lattice spacings. 

The low energy constant g is defined in the chiral 
limit. Therefore, the computations shown in figure [7] 
are extrapolated to m 2 = 0. The final uncertainties are 
dominated by the systematic uncertainty in this step. 
Ref. Il92ll obtained 


g = 0.492(29). 


(115) 


in the chiral limit. 


"Usually, the second term is written as Ofp 2 /(8;r 2 / 2 )) and the ex¬ 
pansion coefficients are assumed to be order one in that variable. With 
a pion decay constant of / ~ Aqcd/4 for our previous estimate of 
Aqcd ~ 500 MeV, there is a numerical difference. Similar ambigui¬ 
ties are present in a prefactor of m up /AQCD, but these numerical fac¬ 
tors do not influence the structure of the expansion and are not well 
known/defined anyway. 


5.2. HQET parameters 

For Nf = 0 the HQET parameters oj, have been de¬ 
termined in |[99l . Here we review the refined strategy 
applied with two flavors of dynamical fermions |[49l . It 
covered the three parameters in the action and the two 
parameters of ,4 | l ) l<-,h ' 1 needed for the computation of / B . 
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Figure 7: A summary of unquenched lattice QCD results for g. The results from the three lattice spacings used in (92) are labeled 'ALPHA'. 
Additionally there are results of Ohki et al. ED. Becirevic et al. El and Detmold et al. El. For Ref. (95), which employs Nf = 2 + 1 dynamical 
flavors, we take the results for a single level of link smearing in the static action. Graph from E3- 


i (t>p CD c.f. large volume limit dominant sensitivity to 



All five used matching observables are similar to what 
we discussed in section [4] We list them in tabled Note 
how three of them are small volume versions of phys¬ 
ical quantities which one would like to predict in large 
volume. If finite volume effects are truly small, these 
quantities are then predicted correctly with very small 
truncation errors of the l/m h expansion. Even with sig¬ 
nificant finite volume effects, this property is expected 
to reduce truncation errors as compared to many other 
choices. 

In static order, only <X>i, <f >2 are needed and they are 
finite volume generalizations of mass and decay con¬ 
stant of the B-meson. In figure [8] we illustrate their de¬ 
termination. The observables are computed for various 
fixed z = ML\ and fixed L\, but different resolutions 
Li/a = 20... 40. They are then extrapolated to the con¬ 
tinuum limit using the asymptotic dominance of ( a/L\ ) 2 
corrections. In order to carry this step out, one needs to 


control the non-perturbative determination of the RGI 
mass M and one needs to know what fixed L means in 
terms of the bare parameters. Indeed, fixed L is only 
defined up to discretisation errors. Here it was cho¬ 
sen to mean fixed gsp(L), the running coupling in the 
Schrodinger functional renormalization scheme. The 
required full control over the renormalization of QCD 
in the light sector had been gained in a series of earlier 
works. Some milestones are l!7^ll00ll771l76l[T0Tlll02l . 
reviewed in some detail in ED- 

The middle column of figure [8] shows the mass- 
dependence of <J), for L = L\ where we match. The 
expected dq = O (M) = O(’) and dH = 0(1) is clearly 
visible, also the 0(1 /M) = 0(1 /z) corrections are seen 
in dH. The right column finally demonstrates the good 
control over the continuum limit in the scaling step to 
Lo, eq. ( |8~i~| i. It is only in this last step where the restric¬ 
tion to the static approximation is relevant at all. Still, 
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the graph does not look much different when the full 
system with five <D, is treated including the NLO cor¬ 
rections ||49| . 

Let us now illustrate these NLO corrections. Figure[9] 
shows the other three observables after subtracting their 
static part. They clearly exhibit the 0( 1 /M) - Of 1 /z) 
behavior at large mass. Again this represents a confir¬ 
mation of the correctness of HQET. 

With these steps carried out, it is only left to evaluate 
eq. (82 1 at the desired lattice spacings where large vol¬ 
ume simulations are carried out. The HQET parameters 
are then known as a function of the RGI mass M, more 
precisely z = ML\ and at a few values of the lattice 
spacing corresponding to integer Li/a. By an interpo¬ 
lation the knowledge of the parameters is extended to 
all values of the lattice spacing within the range corre¬ 
sponding to the accessible Lo/a. The correct value Mb 
for a b-quark in Nature still needs to be determined by 
matching the meson mass as a function of M to the ex¬ 
perimental B-meson mass. This step yields also the first 
prediction of the theory, namely the physical value of 
the renormalized b-quark mass, which is one of the fun¬ 
damental parameters of QCD. We turn to it now. 


5.3. Mass of the b-quark 
5.3.1. Static order 

Lattice HQET computations of the b-quark mass have 
been carried out originally with a just perturbative sub¬ 
traction of the 1 ]a divergence ifTHhl fTQ7l [T()8l in the 
static approximation. A continuum limit does not exist 
in this case, so a corresponding extrapolation may not 
be performed. Instead one can check for stability under 
changes of a. When present, such a stability indicates 
that divergent terms as well as discretisation effects are 
small. For comparison we include a world average of 
the year 2001 by S. Ryan ll 1031 in our summary, table [4] 

The first computations with non-perturbative renor¬ 
malization were still restricted to the quenched approx¬ 
imation. : The initial static order result of O0l was later 
extended to a full NLO mass computation in 11041 . 


5.3.2. HQET at NLO 

Fairly recently a NLO HQET computation with Nf = 
2 flavors of dynamical fermions based on the parame¬ 
ter determination discussed in the previous section was 
completed fl05t . It used the GEVP method {67) for bet¬ 
ter control of the computation of the B-meson mass (at 
a fixed lattice spacing and other bare parameters such 
as //zj’ a,c ) and exploited several CLS lattices. They cover 
three lattice spacings between a = 0.048fm and a = 
0.075fm and several light-quark masses corresponding 


to pion masses 190 MeV < m„ < 440 MeV. At a fixed 
mass of the b-quark, parameterized by z = ML, these re¬ 
sults needed to be extrapolated to the physical pion mass 
and to the continuum limit (all volumes are large enough 
to safely neglect finite volume corrections). This extrap¬ 
olation was performed in the form of one global fit: 

- The non-analytic term in the m 2 expansion. 


is subtracted from the data. 

- The remaining mass-dependence is parameterized 
by a linear term in m 2 n (oc £). 

- The lattice spacing dependence is parameterized 
by the leading term in the Symanzik expansion, 

2 

oc a . 

- The extrapolation is carried out simultaneously for 
two discretizations of HQET (HYP1/2). All fit- 
parameters except for the coefficient of a 2 are com¬ 
mon to both discretizations. 


eq. (|1 11 [), discussed in and known from section 5.1 


- Physical units are taken from a previous extrapo¬ 
lation of the kaon decay constant to the physical 
point and continuum limit ll79l . 


This parameterization of lattice spacing and quark mass 
dependence neglects higher order terms in a 1 and f, also 
the mixed term a 2 x Q It fits the data very well. Data 
and fit, together with the continuum limit, are shown in 
figure[l()]for three prescribed values of z,. 

The results at a = 0 and the physical pion mass were 
then interpolated in z and finally z.b = M h L\ was deter¬ 
mined by requiring 

mnizQ = m“ pcn,ncntal . (116) 


for the interpolation function m B (z), see the right graph 
in figure [T0| The resulting number Mb = 6.58(17)GeV 
is included in our summary, table[4] Repeating the anal¬ 
ysis with all 1//Hh terms dropped yields the static result. 
It is almost identical indicating that NNLO terms are 
very small. Note that a naive estimate of these terms, 
A 3 /m^ ~ 4MeV, is indeed very small compared to the 
dominating statistical errors. 

The RGI mass Mb can be changed easily to /«* de¬ 
fined earlier, using a knowledge of the A-parameter 
ED and inserting the MS 4-loop r and f3 functions 
EOEUHIEDI int0 eq- © and eq. |p). 

These numbers (last row of table^f are well con¬ 
trolled. Given the small difference to Nf = 0 results, 
it is even very hard to imagine that the missing strange 
quark is a significant source of error. The numbers agree 
with other determinations and the PDG average II1091 of 
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Figure 8: Continuum extrapolation of Oi and <&2 in L = L\ (left). Their resulting mass-dependence (z = ML, middle). On the right we show the 
step scaling to obtain them in L = L 2 in LO HQET. The b-quark in Nature has about 1/z « 0.08. 


/n b (w b )[GeV], MS M b [GeV] 


N f 

ref. 

remarks 

static 

0(1 /uih) 

static 

0(1/m h ) 

0-2 

H03| 

2-loop subtracted; 2001 average 

4.30(10) 




0 

El 

NP 

4.12(7)(4) 




0 

|104| 

NP 

4.32( 5) 

4.35( 5) 

6.81( 8) 

6.76( 9) 

2 

|105| 

NP 

4.21(11) 

4.21(11) 

6.57(17) 

6.58(17) 


Table 4: HQET results for the RGI b-quark mass Mb and m* = rab(mb) in the MS scheme computed using lattice HQET. Results with perturbative 
subtraction of 1 /a divergences are labeled “n-loop subtracted”. When renormalization and matching is treated non-perturbatively, we just have a 
label “NP”. 


m b * = 4.18(3). They thus boost our confidence in hav¬ 
ing systematic errors in the determination of m b under 
control. 


5.4. B-meson decay constants 
5.4.1. Mass-dependence 

Before discussing the most up-to-date Nf = 2 results 
for the decay constants, we spend some time on an im¬ 
portant theoretical aspect: what is the dependence of 
the decay constants on the heavy quark mass? Does 
the asymptotic mass-scaling eq. ( |7T| even reach down 
to m h = m c and how far is the b-quark from the heavy 


quark limit? With good statistical and systematic pre¬ 
cision, these questions have only been studied in the 
quenched approximation [ 1041 ; 110 1. For qualitative ear¬ 
lier results we refer to l fTl2 1031 . In 11101 the static 
matrix elements were first determined at three differ¬ 
ent lattice spacings increasing the precision by a use 
of translation invariance as well as the GEVP method. 
They were then extrapolated to the continuum limit as 
shown in figure 11 Only the B s meson was investi¬ 
gated because the chiral limit is singular in the quenched 
approximation, rendering computations of the limit of 
very light quarks rather meaningless. 


The static result for the ground state is compared to 
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Figure 9: oj 1 /,Wh) (Li) = Op CD (Li) - i//(Li), / = 3,4 in L\ on the left and their step scaling function to obtain = <I> 9 CD (L 2 ) - 77 /(L 2 ) on 

the right. The bottom line shows ® 5 (Li); spin symmetry means that its static part vanishes, 775 = 0. The b-quark in Nature has about l/z « 0.08. 
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Figure 10: Left: Chiral and continuum extrapolation of mg for the z -values used in the determination of Zb■ Open/filled symbols refer to 
HYP1/HYP2 data points as do long/short dashed curves, respectively. The continuum pion mass dependence is given by the solid lines together 
with a shaded error band. The triangle shows the estimated continuum limit at the physical pion mass. The figure uses a notation y = £ and 6 is an 
index for the different disretizations. 

Right: Interpolation to Zb by imposing eq. jl 16) . Graphs from tm 





/ B [MeV] 

/b s [MeV] 

hJh 

Nf 

ref. 

remarks 

static 0(l/wh) 

static 0(l/«Jh) 

static 0( 1 / mh) 

0 

11121 

summary 

PT 

276(55)(19) 


1.22(4) (2) 

0 

[110| 

NP 


229(3) 216(5) 


2 

Dual 

NP 

190(5)(2) 186(13) 

226(6)(9) 224(14) 

1.189(24)(30) 1.203(65) 

2+1 

1:1141 

PT, 

1-loop 

219(17) 

264(19) 

1.193(41) 

2 

im 

FLAG 

189(8) 

228(8) 

1.206(24) 



average 1 




2+1 

tna 

FLAG 

190.5(4.2) 

227.7(4.5) 

1.202(22) 



average 1 





Table 5: Results for the B-meson decay constant computed using lattice HQET in the upper part. Results with with completely perturbative 
renormalization are labeled by “PT”. When renormalization and matching is treated non-perturbatively, we just label “NP”. 

1 The FLAG averages in the lower part combine various methods, not just HQET. 
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fi(y,a )/GeV fi(y,a )/GeV 



Figure 13: Extrapolation of the Z? (left panel) and Z? s (right panel) meson decay constant to the physical point. On the left, the extrapolation using 
HM^-PT at NLO (filled triangle) is compared to a linear one (open triangle), in order to extract the systematic error from truncating HM^PT at 
NLO. For /b s only a LO formula is known and shown. As a comparison also the final result, the continuum value of /b s = [/b s //b]/b is shown. 
Figure from nm 



Figure 11: Continuum limit of the Nf = 0 static matrix element = 
(0| AQ tat |fi) in RGI normalization. The lower part shows the ground 
state matrix element and the upper part the first excited state with the 
quantum number of the B -meson. Graph from firm 



1/( r o m PS ) 


Figure 12: A^f = 0 heavy quark mass dependence of the heavy-strange 
decay constant. All dimension-full quantities are expressed in units of 
r 0 a 0.5fm fiTO . The charm quark is located at about rapsro = 0.2 
and the b-quark at mps/"o = 0.07. Graph from htqi . 


results at finite mass in figure 12 The graph shows 
/ps y/mps/Cps, where PS is a strange-heavy pseudo¬ 
scalar against the inverse heavy meson mass. At the 
lowest order, the meson mass is proportional to the 
quark mass. Therefore mps can be taken as a proxy for 
the quark mass. The conversion function Cps is com¬ 
puted with the full perturbative knowledge, i.e. it has 
relative errors of order cr(/«h) 4 - However, as discussed 
before, for b- or c- quark masses it is not at all clear that 
the perturbative errors are negligible. 


The graph shows a very consistent picture, with the 


static limit computation on the left being in line with the 
relativistic data points on the right and an HQET com¬ 
putation with non-perturbative matching and including 
the 1 /mh terms at mp$ = nip. Given the somewhat un¬ 
certain status of the perturbative Cps, the good agree¬ 
ment with a simple linear behavior in 1/nips ~ l/m^ 
over a large range is in fact somewhat surprising. 

In the above discussion of the mass-dependence it 
is hard to get around a perturbatively extracted CpsPl 


'-Non-perturbative definitions of Cps do not separate the logarith- 
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Without it taking care of the logarithmic modifications 
of the naive scaling law the curve in figure [12] would 
not even have a finite limit at 1/wh = 0. However, 
for predictions at the b-quark mass (or any other large 
mass), the non-perturbative HQET parameters w, from 
section 5.2 can be used. From now on we only discuss 


ALPHA collaboration results obtained in this entirely 
non-perturbative way. Only in the comparison to results 
in the literature perturbatively renormalized results by 
other groups are shown. 


5.4.2. HQET results with non-perturbative parameters 
The physical value of the b-quark mass (parameter¬ 
ized by Zb) is known from the previous section. All 
HQET parameters are then interpolated to that quark 
mass. Once they are known, the computation of the de¬ 
cay constants is in principle straight forward. However, 
again care has to be taken about several limits. 


■ to -2 



Figure 14: Plateau averages after applying the GEVP analysis to 
data obtained on the Nf = 2 CLS ensemble N6 (a = 0.048 fm, 
m n = 340 MeV). The upper plot shows the B s -meson static matrix 
element p s tat- The lower plot shows the chromo-magnetic matrix ele¬ 
ment p S pi n . Note eq. jl 17k for how these terms enter in the expression 
for the decay constant. Figure from urn 


1) for each lattice spacing, volume and light quark 
mass, the decay constant has to be extracted from the 
large-time limit of correlation functions. One usually 
forms ratios of correlation functions which have the de¬ 
cay constant as their large time (= t ) limit and asymp¬ 
totic corrections to it of ordcipj A/A.i t exp(-A/A,i t). 


mic corrections from power corrections. 

13 The pre-factor A E t is present since we are dealing with an effec¬ 
tive theory computation, see rrm 


Here, t is the time extent of the two-point functions en¬ 
tering and A E„i — E n - E\ is the difference of the in¬ 
dicated energy levels; AE^.i is around 600 MeV. Alter¬ 
natively, with the GEVP method, the dominant correc¬ 
tions can be changed to 0(AE, h i to exp(-A E n \ to)) with 
a larger n. At the same time, a second time separation 
to is present in the GEVP which in practice is to = f/2. 
For a proper explanation we have to refer to ll67l 11131 . 
The analysis in the latter reference estimates the correc¬ 
tion term from the results at smaller f (l and then uses to 
large enough such that the estimated correction is a fac¬ 
tor three below the statistical error of the final result. We 
show two such plateau analysis for the B s meson in fig¬ 
ure 14 The plotted quantities p spm enter the decay 


constant through the HQET expansion 


/b, = exp (x)/^a 3 m B J2, 

X = ln(Z^ QET ) + ln(a 3/2 pstat) (117) 

Pkin ~t" tU S pi n /?spin O/\ 0l p^ (l) . 

The chosen start of the plateau averages according to the 
explained criterion look overly conservative in the plot, 
but on the other hand, excited state contaminations are 
clearly visible in the static piece, /? stat , and they are also 
hard to rule out for the 1 /ni tl term shown. 

In summary, obtaining the ground state matrix ele¬ 
ment is far from trivial, but within the rather conserva¬ 
tive errors of 111 13l it appears to be under control. We 
do not want to hide that the statistical errors for the B- 
meson (rather than B s ) are somewhat larger, see Figure 

i of nm 

2) A continuum extrapolation has to be carried out 
and 

3) the results have to be extrapolated to the physical 
quark mass. 

Just like in the case of the meson mass, 2) and 3) are 
carried out in one global fit. As seen in figure [13] the 
lattice spacing dependence is not significant at all and 
one could just average the numbers at the different lat¬ 
tice spacings. However, this would not account for the 
uncertainty in the statement that the a-dependence is in¬ 
significant. Therefore, as above, an a 2 term is fitted. The 
light-quark-mass dependence is parameterized as a lin¬ 
ear dependence and alternatively as the one predicted by 
heavy meson chiral perturbation theory. The difference 
in the extrapolated value is very small and accounted for 
in the errors. 

4) In principle we should also worry about effects of 
the finite volume, but these are very small for the con¬ 
sidered volumes. 

The resulting numbers are listed in table [5] together 
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with the previously computed Nf = 0 decay con¬ 
stants.^] The latter were quoted as tm 

= 2.14(4), (118) 

r l /2 /£ QET V^ = 2.02(5). (119) 

We combine them with the present estimate ro = 
0.50(2) fm (see JTT7]) to fgf = 229(3)(10) MeV and 
^hqet _ 216(5)(9) MeV. In the table we also show 
older numbers as well as the recent computation by the 
RBC/UKQCD collaboration [TT1Q3H- AH of these 
use a perturbative renormalization. The apparent dif¬ 
ference between 2+1 and 2 flavor results in the static 
approximation are likely due to the different renormal¬ 
ization and matching. We remind the reader that static 
results have a dependence on the matching condition 
which is of order 1/m h. It would therefore be prema¬ 
ture to conclude that the difference between the static 
numbers in the table is due to perturbative renormaliza¬ 
tion. 

As a separate result, it was also observed, in the con¬ 
tinuum limit of the Nf - 0 theory, that the decay con¬ 
stant of the excited state is larger than the ground state 
decay constant by the following amount II 101 : 


rstat 
^ B' 


f stat 
■'Bs 



1.24(7) for Nf - 0. 


( 120 ) 


A careful continuum limit was needed to observe this 
splitting of decay constants, see figure [XT] Still it is in 
agreement with an older investigation where a-effects 
were not yet controlled Ml 1911 . 


6. Conclusions, opportunities and challenges 

So far, numerical results of non-perturbative HQET 
have been obtained only in the theory with at most 
two dynamical quark flavors. The dominant impact 
of the project described here concerns the concept, the 
methodology and the qualitative features of the HQET 
expansion. It has been shown that this works out, in 
principle and in practice. Having said that, we would 
like to emphasize, however, that the numbers in table [4] 
and table [5] provide a very valuable crosscheck on other 
results for flavor physics, even if they are obtained with 
just up and down quarks as dynamical quarks. The rea¬ 
son is that there is overwhelming numerical evidence 


14 Earlier Nf = 0 static results are summarised in ED:- A later 
static decay constant at a single lattice spacing is found in fTin . The 
results of fnol are much more precise and do have a continuum limit 
extrapolation. 


3.0 

2.5 

2.0 

c 

< 1.5 
1.0 
0.5 
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f4 
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— HYP1 
HYP2 


0.000 0.002 0.004 0.006 0.008 

a 2 [fm 2 ] 


Figure 15: Continuum limits (stars) of AE^ (circles) and AZsi^ 1 ^ 
(diamonds). Shown are the results for both HYP1 (red, shifted to the 
left) and HYP2 (blue, shifted to the right). Graph from 11211 . 


that strange quark loop effects are strongly suppressed 
when one looks at the low energy properties of QCD 
mam Therefore, it is very reasonable to expect 
present uncertainties to be dominated by the errors in 
renormalization and matching or in the assumptions 
made on the size of lattice spacing effects or the associ¬ 
ated extrapolations. In these respects, HQET is on very 
solid theoretical grounds. There are no doubts about 
its renormalizability and the existence of the continuum 
limit, even if these properties are not proven to all orders 
of perturbations theory (see section |1.6| >. Renormal¬ 
ization and matching is carried out non-perturbatively 
without any compromise. 

A crucial question concerning any expansion is of 
course how large truncation errors are. First of all, we 
do not expect the HQET expansion to be a convergent 
one. This does not matter, good asymptotic expansions 
are just fine, especially if one is able to compute only a 
few terms anyway. Looking at the size of the computed 
corrections, it is important to remember that they do de¬ 
pend on the matching conditions (see section |23) . The 
way we (mostly) perform the matching, HQET is used 
to compute the finite size effects between e.g. a finite 
volume definition of a decay constant and the true large 
volume one. Apart from the fact that for B-physics we 
are expanding in a truly small parameter A/m\, ~ 1/10, 
this is one explanation for the very small NLO correc¬ 
tions that were found. However, there are cases where 
a priory we do not have a reason to expect such addi¬ 
tional suppressions due to the matching condition. One 
such case is the splitting between ground state and first 
(“radially”) excited state of a B-meson 11211 . Its con- 
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tinuum limit in the quenched approximation is shown in 
the lower part of figure [15] It is remarkable that there is 
no significant difference between the pure static result 
and the one including 1 /m h corrections after the contin¬ 
uum limit has been taken. In numbers, the static split¬ 
ting is r 0 A£ff = 1.50(5) while r 0 AE ^ /m) = 0.03(6) 
which is a bit smaller than expected from a simple 
1.50 x A/mb ~ 0.15 estimate. 

There are good opportunities for improving the 
present numerical results in the near future. They arise 
because the uncertainties are (apart from the truncation 
error which we discussed is small) by far dominated 
by statistical errors apart from, maybe, the one coming 
from the omission of the strange quark vacuum polar¬ 
ization. The latter will be removed by the new set of 
large volume CLS simulations ll 1221 (and the associated 
matching program) and there is also plenty of room to 
enhance the precision in the finite volume matching and 
step scaling. 

There is also a very good opportunity to help in the 
understanding of the so-called V U b puzzle, which says 
that this CKM matrix element differs between different 
determinations. Form factors of semi-leptonic B-decays 
are a key to understanding the puzzle. Work on them 
requires an extended program to match the full set of 
currents BTl [52] [53l [54) as well as the computation of 
form factors in large volume HQET I123II1241I125I . We 
refer to [55) for a review. 


Let us finish with a challenge for the future. It is 
the exponential deterioration of the signal-to-noise ra¬ 
tio, eq. ©. There are two issues. The first is that E stat 
has a linear divergence which means the problem be¬ 
comes worse as we approach the continuum limit fur¬ 
ther. As we have explained in section[3]this issue is not 
so severe since the coefficient of the 1/a divergence is 
very small for the HYP2 static action. However, the fi¬ 
nite part is not small and leads to rather short plateaus, 
see figure |~[4| Formula (91 1 holds for our standard esti¬ 
mators of the correlation functions and the standard im¬ 
portance sampling. We should not give up to think about 
whether there are ways around it and see long plateaus 
with small errors in the future. 
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